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ABSTRACT 

We investigate the gravitational wake due to, and dynamical friction on, a 
perturber moving on a circular orbit in a uniform gaseous medium using a semi- 
analytic method. This work is a straightforward extension of Ostriker (1999) 
who studied the case of a straight-line trajectory. The circular orbit causes 
the bending of the wake in the background medium along the orbit, forming a 
long trailing tail. The wake distribution is thus asymmetric, giving rise to the 
drag forces in both opposite (azimuthal) and lateral (radial) directions to the 
motion of the perturber, although the latter does not contribute to orbital decay 
much. For subsonic motion, the density wake with a weak tail is simply a curved 
version of that in Ostriker and does not exhibit the front-back symmetry. The 
resulting drag force in the opposite direction is remarkably similar to the finite- 
time, linear-trajectory counterpart. On the other hand, a supersonic perturber 
is able to overtake its own wake, possibly multiple times, and develops a very 
pronounced tail. The supersonic tail surrounds the perturber in a trailing spiral 
fashion, enhancing the perturbed density at the back as well as far front of the 
perturber. We provide the fitting formulae for the drag forces as functions of the 
Mach number, whose azimuthal part is surprisingly in good agreement with the 
Ostriker's formula, provided Vpt = 2Rp, where Vp and Rp are the velocity and 
orbital radius of the perturber, respectively. 

Subject headings: hydrodynamics — galaxies : kinematics and dynamics — ISM: 
general — shock waves 

1. INTRODUCTION 

Dynamical friction refers to momentum loss suffered by a massive object moving through 
a background medium due to its gravitational interaction with its own induced wake. The 
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gravitational drag force removes angular momentum from an object in orbital motion, caus- 
ing it to gradually spiral in toward the center of the orbit. In a pioneering study, Chan- 
drasekhar (1943) derived the classical formula of dynamical friction in a uniform collisionless 
background, the result of which has been applied to a number of astronomical systems. Ex- 
amples include orbital decay of satellite galaxies orbiting their host galaxies (e.g., Tremaine 
1976; Lin & Tremaine 1983; Weinberg 1989; Hashimoto et al. 2003; Fujii et al. 2006; see 
also Binney & Tremaine 1987), dynamical fates of globular clusters near the Galactic cen- 
ter (e.g., Kim & Morris 2003; McMillan & Portegies Zwart 2003; Kim et al. 2004), galaxy 
formation within the framework of hierarchical clustering scenario (e.g., Zentner & Bullock 
2003; Bullock & Johnston 2005 and references therein), formation of Kuiper-belt binaries 
(Goldrcich ct al. 2002), and planet migrations (Del Popolo et al. 2003) via interactions with 
planetesimals, etc. 

Although less well recognized, dynamical friction also operates in gaseous backgrounds. 
Using a time-dependent linear perturbation theory, Ostriker (1999, hereafter 099) derived 
the analytic expressions for the density wake and drag force for a perturbcr in a uniform 
gaseous medium. 099 showed that resonant interactions between a perturber and pressure 
waves make the gaseous drag more efficient than the collisionless drag when the Mach number 
A4 ^ 1. She also found that even a subsonic perturber experiences a nonvanishing gaseous 
drag if interaction time between the perturber and the background is finite. This is an 
improvement on the previous notion that the gaseous drag is absent for subsonic perturbers 
because of the front-back symmetry in the steady-state density wake (Dokuchaev 1964; 
Ruderman & Spiegel 1971; Rephaeli & Salpeter 1980). 

The results of 099 were confirmed numerically by Sanchez-Salcedo & Brandenburg 
(1999) and have been applied to various situations including massive black hole mergers 
in galactic nuclei (e.g., Escala et al. 2004; Liu 2004; Escala et al. 2005; Dotti et al. 2006), 
orbital decay of compact objects (e.g., Narayan 2000; Karas & Subr 2001) and associated 
viscous heating (e.g., Chang 2001; Chang & Choi 2003) in accretion disks, and heating of 
an intraclustcr medium by supersonically moving galaxies in clusters (e.g., El-Zant ct al. 
2004; Faltenbacher et al. 2005; Kim et al. 2005). Without involving shocks, density wakes 
in a collisionless medium are distributed more smoothly and achieve larger amplitudes than 
those in a gaseous medium (Mulder 1983), which led Furlanetto & Loeb (2002) to suggest 
that X-ray emissions from galaxy wakes can in principle be used to discern the coUisional 
character of dark matter in galaxy clusters. 

While the results of Chandrasekhar (1943) and 099 are simple and provide good physical 
insights, they apply strictly to a mass traveling on a straight-line trajectory through an 
infinite homogeneous background. Real astronomical systems obviously have nonuniform 
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density distributions and perturbers tend to follow curvilinear orbits. For instance, motions 
of galaxies in galaxy clusters, binary black holes near the central parts of galaxies, and 
compact stars in accretion disks are better approximated by near-circular than straight-line 
orbits, and their background media usually in hydrostatic equilibrium are stratified in the 
radial direction. Even for objects experiencing orbital decay, a near-circular orbit is a good 
approximation if the associated friction time is longer than the orbital time. Consideration 
of a circular-orbit perturber is of particular interest since it will allow the perturber, if 
supersonic, to overtake the backside of its wake that was created about an orbital period 
earlier. In this steady-state wake that eventually forms has morphology and drag 

force that might be significantly different from the linear-trajectory counterparts. Numerical 
simulations carried out by Sanchez-Salcedo & Brandenburg (2001) and Escala et al. (2004) 
indeed show that the density wake by a near-circular orbit perturber contains a trailing spiral 
tail, which is absent in the linear-trajectory cases. They also showed that the resulting drag 
force is smaller than the estimate based on the formula given in 099, which is probably due 
to the near-circular orbit, although the effect of nonuniform backgrounds in their models 
cannot be ignored completely. 

The drag formula based on perturbers moving straight in either a collisionless medium 
or a gaseous medium depends on the Coulomb logarithm ln(A) = ln(rmax/Tmin), where rmin 
and Tmax are the cutoff radii introduced to avoid a divergence of the force integrals.^ While 
many previous studies conventionally adopted rmin and r^ax as the characteristic sizes of the 
perturber and the background medium, respectively, the choice of r^ax remains somewhat 
ambiguous for objects moving on near-circular orbits (e.g., Binney & Tremaine 1987). For 
collisionless backgrounds, Hashimoto et al. (2003) and Fujii et al. (2006) performed N-body 
experiments for orbital evolution of satellite galaxies in a spherical halo, and found that 
the Coulomb logarithm with rmax varying proportionally to the orbital radius rather than 
fixed to the system size gives better fits to their numerical results (see also Tremaine 1976). 
We shall show in the present work that a similar modification of the Coulomb logarithm is 
necessary in order to apply the results of 099 to perturbers on near-circular orbits in gaseous 
backgrounds, as well. 

In this paper, we consider a perturber moving on a circular orbit in a uniform gaseous 
medium. Using a linear semi-analytic approach, we explore the structure of the density wake, 
evaluate the drag force on the perturber, and compare them with those in the straight-line 
trajectory cases. In §2, we revisit the linear perturbation analysis of 099 for the perturbed 
density response and apply it to the case of a circular-orbit perturber. We solve the resulting 



For a perturber moving with velocity Vp through a gaseous medium, rmax = Vpt, where t denotes time 
elapsed since the perturber was introduced (099). 
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equations numerically. In §3, we present the numerical results for the wakes and drag forces 
with varying Mach number. We provide simple fitting expressions to the numerical results, 
and show that ©striker's formula still gives a good estimate for the drag force on a circular- 
orbit perturber only if the outer cutoff radius in the Coulomb logarithm is taken equal to the 
orbital diameter of the perturber. In §4, we summarize the present work and briefly discuss 
our finding. 



2. FORMULATION 

2.1. Formal Solution for Density Wake 

We consider the response of gas to a point-mass perturber moving on a circular orbit 
and calculate the resulting gravitational drag force on the perturber. We treat the gas as an 
inviscid, adiabatic fluid, and do not consider the effects of magnetic fields as well as gaseous 
self-gravity. The governing equations for ideal hydrodynamics are 

|^ + V-(pt;) = 0, (1) 

and 

^ + VVV = —VP - V$ext, (2) 

ot p 

where $ext is the gravitational potential of the perturber. Other symbols have their usual 
meanings. 

Following 099, we consider an initially uniform gaseous medium with density po- As- 
suming that the wake induced by the perturber remains at small amplitudes, we linearize 
equations (1) and (2) using p = Pq[1 + a{x,t)] and v = Cs/3{x,t), where Cg is the adiabatic 
speed of sound in the unperturbed medium and a and f3 denote the dimensionless density 
and velocity perturbations, respectively. Eliminating (3 from the linearized equations, one 
obtains a three-dimensional wave equation 



^ « - = -^Pext(a;,t), (3) 



where pext = V^$ext/(47rG) represents the mass density of the perturber. The formal solu- 
tion to equation (3) based on the the retarded Green function technique is given by 

a(x, t)^-JI d'x'dt' pext(a^', (4) 

Cc J J \x \ 
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(099; see also Jackson 1975). In the case of a perturber on a straight-line trajectory, 099 
solved equation (4) directly to obtain an expression for the perturbed density. The same 
result was found by Furlanetto & Loeb (2002) who independently used a Fourier transform 
method in both space and time variables. 



2.2. Density Wake for Circular-Orbit Perturbers 

We now concentrate on the case where a point-mass perturber with mass Mp moves on 
a circular orbit with a fixed orbital radius Rp and a constant velocity Vp in an otherwise 
uniform gaseous medium; the angular speed of the perturber is Q = Vp/ Rp. It is convenient 
to work in cylindrical coordinates (i?, (/?, z) whose origin lies at the center of the orbit. The 
z-axis points perpendicular to the orbital plane. Assuming that the perturber is introduced 
at (i?p,0,0) when t = 0, one can write pext(ic,t) = Mp 6{R — Rp)6[Rp{ip — Qt)]6{z)T-C{t), 
where H{t) is a Heaviside step function. Equation (4) is then reduced to 

ciRp J d{w]R,z) \ n J 

where w = tp' — tp and s = ip — VLt are angular distances in the z = Q plane^, 

d{w;R,z)=^- — - — ^= l + i?^ + f^-2i?cosw , (6) 
Rp \ / 

and M. = Vp/cs is the Mach number of the perturber. In equation (6), R = R/Rp and 
z = z/Rp. 

Figure 1 schematically illustrates the situation at the orbital plane and the meanings 
of variables used in equation (5). At time the perturber is located at Xp — {Rp,^t,0). 
During its journey along the thick curve, the perturber continuously launches sound waves 
that propagate into the background gaseous medium. The position x = {R, ip, z) denotes a 
region of interest in the surrounding gas where the density response will be evaluated. Since 
the sound waves have finite traveling time, only the signals emitted by the perturber at the 
location(s) x' = {Rp,ip',0) at the retarded time f = t — \x — x'\/cs are able to affect the 
point x at time t. Note that s and w represent the projected angular distances in the orbital 
plane between Xp and x and between x' and x, respectively. The symbol d in equation (6) 
refers to the three-dimensional linear distance between x' and x normalized by Rp. 



^Note that w and s in 099 are defined as linear distances along the line of motion, while they measure 
angular distances in the present work. 
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Using the identity 5 [f{w)] — S{w—Wi) /\f'{wi) \ , where Wi are the roots of an arbitrary 
function f{w), equation (5) is further simphfied to 

a{x,t)^^V{x,t), (7) 

with the dimensionless perturbed density 

vix, = V n ( ^] . (8) 

\wi + s - M'^R sinwA \ il J 
Here, the summation is over all possible roots Wi that satisfy the condition 

M d{wi;R,z) = -{wi + s), (9) 

for fixed values of R, z, and s. The function T-C[{wi + <f)/^] in equation (8) defines the 
region of influence (or casual region) outside of which sonic perturbations sent off by the 
perturber at t = have insufficient time to reach. For R,z,flt 3> 1, equation (9) yields 
Wi + (p — Qt — A4{R^ + z'^y^'^, so that the region of influence corresponds roughly to a sphere 
with radius Cst centered at the orbital center. Appendix A presents limiting solutions of 
equation (8) near the perturber. 

Since o? is a periodic function of w with period 27r, equation (9) has at least one real 
root and may possess multiple roots for Wi depending on the values of A^, R, z, and s. In 
Appendix B, we describe how the number of solutions of equation (9) vary with the Mach 
number of a perturber. It turns out that there is only a single root for Wi everywhere for a 
subsonic, circular-orbit perturber, which is the same as in the straight-line trajectory case 
(099). When a circular-orbit perturber moves at a supersonic speed, however, equation (9) 
has an odd number of roots that contribute to the wake in a steady state, which is distinct 
from the straight-line trajectory case where only one or two points along the orbit influence 
the wake. 



2.3. Gravitational Drag Force 

Once the gravitational wake a{x,t) is found, it is straightforward to evaluate the drag 
force exerted on the perturber: 

Fdf = GM.po [ d'x »M^^-^p) _ (10) 

J I X Xp I 

In the straight-line trajectory case studied by 099, a{x, t) always remains cylindrically 
symmetric with respect to the line of motion, resulting in the drag force in the anti-parallel 
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direction. When the perturber is on a circular orbit, a{x,t) loses the cylindrical symmetry 
and instead becomes symmetric relative to the orbital plane, making the vertical component 
of i^DF vanish. 

We decompose the nonvanishing parts into the radial and azimuthal components: 

For = -^(XRii + I.*). ^.i^^^H^ (11) 

p 

where 
and 

I ^-Ml f rf3^ Vix,t)Rsins 

47r 7 (l + F + i?2_2i?coss)3/2 ^ ^ 

Note that measures the drag force along the lateral direction of the instantaneous per- 
turber motion, while X^ is for the backward direction. The dimensional term T in equation 
(11) allows to directly compare Xr andX^ with the linear-trajectory counterparts (see eq. [12] 
of 099). As we shall show in §4, it is the azimuthal drag X^ that is responsible for the orbital 
decay of a perturber. 



2.4. Numerical Method 

We solve equations (8) and (9) numerically to find the perturbed density distribution 
T>{x, t) for given A4 and t. We first construct a three-dimensional Cartesian mesh centered at 
the center of the orbit, and solve equation (9) for Wi at each grid point using a hybrid Newton- 
bisection method. By checking the conditions for multiple roots discussed in Appendix B, 
we ensure that we do not miss any solution for Wi. The corresponding drag forces Xr and 
X^ are calculated by direct integration of equations (12). 

Since the density wake often exhibits sharp discontinuities especially for supersonic 
perturbers and is distributed over a large spatial range, it is important to check that the drag 
forces we calculate do not depend on the size of the computational box and its resolution. For 
fixed Ai, we repeated the calculations with varying box size and resolution and found that 
depending on the Mach number, the box size of ~ (20 — 100)-Rp and resolution of ~ (80 — 640) 
grids per IRp are sufficient to guarantee good convergence of the drag forces. Although the 
density perturbations are non-zero outside the box, they have very low amplitudes and are 
located far from the perturber, providing a negligible contribution to the drag. Very high 
resolution calculations are required for Mach numbers near the critical values A4n, in which 
cases the wake tails become thin and dense (see Appendix B). 
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3. 



RESULTS 



3.1. 



Density Wake 



3.1.1. 



Supersonic Cases 



We begin by illustrating temporal evolution of density perturbations induced by a super- 
sonic perturber. Figure 2 shows snapshots of the density wake and the corresponding number 
of roots of equation (9) at the orbital plane for Ai = 2.0. Time is expressed in units of Rp/cg. 
The black circle represents the orbit of a perturber which, introduced at (i?, 93, z) = {Rp, 0, 0) 
initially, moves in the counterclockwise direction. At early time {t ^ 1.5), the density wake 
consists of a Mach cone and a sonic sphere that are curved along the orbit. Except the 
bending, the overall wake structure, not to mention the number of solutions for Wi which is 
either one or two inside the casual region, is the same as in the straight-line trajectory case 
of 099. As the wake bends further, the Mach cone and the sonic sphere become folded at 
the innermost interface, creating high-density regions near the center (t = 1.8). The wake 
expands with time at a sonic speed and the Mach cone becomes elongated further. 

Unlike in the case of a straight- line trajectory where the Mach cone and the sonic sphere 
never interact with each other, the perturber (and the head of the Mach cone) on a circular 
orbit is able to enter its own wake, providing additional perturbations for some regions inside 
the sonic sphere. Alternatively, this can be viewed as the sonic sphere whose center lies at the 
initial position of the perturber expands radially outward, swallowing a part of the elongated 
Mach cone. Consequently, the overlapping of the Mach cone and the sonic sphere creates 
a high-density trailing tail that has received perturbations three times from the perturber 
{t ^ 1.9). Figure 2 shows that immediate outside the sonic sphere, there still exists a region 
of the undisturbed Mach cone with two WiS. As time proceeds, however, this region moves 
away from the perturber and thus gives an increasingly small contribution to the drag force. 
In a steady state which is attained at i ^ 00, the entire domain is affected either once or 
three times by the perturber. 

Figure 3 displays the steady-state distributions of the density wake for At = 2.0 on the 
x/Rp = 1, y = 0, and z = planes, which clearly shows that the trailing tail loosely wraps 
around the perturber. The tail becomes narrower with increasing \z\. A close inspection 
of the tail at the z = plane shows that density becomes smaller away from the perturber 
along the tail and is largest at the edges across the tail. This is similar to the density 
distribution inside the Mach cone in the linear trajectory case where diverging flows (and 
reduced gravitational potential) after the shock cause density to decrease away from the 
shock front (and the perturber). This suggests that the edges of the tail are shock fronts. 
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The outer edge of the tail shown in Figure 3 that connects smoothly all the way to the 
perturber corresponds indeed to the surface of the curved Mach cone (see also Fig. 1). On 
the other hand, the inner edge of the tail traces the interface between the Mach cone and the 
sonic sphere that is newly swept up by the expanding sonic sphere (see Fig. 2). Appendix C 
proves that the half-opening angle of the head of the curved Mach cone in the z = plane 
is equal to sin~^(l/A1), entirely consistent with the case of a linear trajectory. 

The wake tail becomes thicker as Ai increases from unity. Figure 4 shows the steady- 
state density wakes as well as the number of -ujj's that contribute to V ioi M. = 4.0 and 
M. = 5.0 at the z = plane. In both panels, the perturber moving in the counterclockwise 
direction is located at x/Rp = l,y = 0. When M. = 4.0, the tail is fat enough to cover most 
of the space except near the center and a narrow lane between the tail edges. At Mi ~ 4.603, 
the inner edge of the tail eventually touches the neighboring outer edge, enabling three Wi 
for the entire region under the conditions expressed by equation (Bl). When the Mach 
number is slightly larger than Aii, the tail overlaps itself. This in turn creates a new tail 
with five wi's. as Figure Ah displays. As A4 increases further, the new tail again becomes 
thicker, starts to overlap at Ai2 ~ 7.790, and produces a narrow lane with seven Wj's when 
> Ai2- The same pattern repeats with increasing J\4, and equation (B3) determines the 
critical Mach numbers. 



3.1.2. Subsonic Cases 

Unlike in the supersonic cases where a perturber can overtake its own wake and create 
a tail with complicated structure, sonic perturbations generated by a subsonic perturber 
produce a gravitational wake that is spatially smooth and does not involve a shock. Since 
perturbations propagate faster than a perturber with Ai < 1, the whole casual region is 
affected by a perturber just once, corresponding to a single Wi at any position. Figure 5 shows 
the slices of the perturbed density for M. — 0.5 in the x/Rp — 1, y — 0, and z — planes 
when a steady state is reached. Again, the perturber is located ai x/ Rp — 1, y — z — 0. The 
density structure in the z — plane is simply a curved version of what a linear-trajectory 
perturber would produce. In particular, as Appendix A shows, the iso-density surfaces near 
the perturber have the shapes of oblate spheroids with ellipticity e = Ai, with the short axis 
parallel to the direction of the motion of the perturber, a characteristic feature of a subsonic 
wake created by a linear-trajectory perturber (099). 

Notice, however, that the bending of wakes in circular-orbit cases makes the perturbed 
density distributions intrinsically asymmetric. This results in nonvanishing drag forces even 
in a steady state, and the dominant contribution to the drag comes from high-density regions 
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near the perturber. This is markedly different from the purely steady-state linear-trajectory 

cases where a subsonic perturber experiences no drag due to the front-back symmetry of a 
wake about the perturber (e.g., Rephaeli & Salpeter 1980). Even if the finite interaction 
time between the straight-line perturber and the background gas is considered, regions with 
symmetric perturbed density close to the perturber exert zero net force (099). Nevertheless, 
the resulting drag force in the backward direction of motion on a circular-orbit perturber is 
almost the same as that in the linear-trajectory cases, as we will show in the next subsection. 
Compared with supersonic cases, the tail in a subsonic wake is short, loosely wound, and 
very weak, suggesting that its contribution to the drag force is negligible. 



3.2. Gravitational Drag Force 

As sonic perturbations launched from a perturber at t = propagate radially outward, 
the volume of space exerting the gravitational drag on the perturber grows with time. Figure 
6 plots the drag force as functions of time for a few chosen Mach numbers. The solid and 
dotted fines are for X^p and Xr, respectively. Since V{x,t) is singular x — Xp, only the 
region with r > rmin = Rp/ 10 is taken into account in force computation, where r is the 
three-dimensional distance measured from the perturber; the dependence on rmin will be 
checked below. The drag force on a subsonic perturber increases almost linearly witfi time 
before turning abruptly to a constant value, whereas a supersonic wake with a high-density 
tail gives rise to slow fiuctuations in tfie drag at early time. At any event, botfi components of 
the drag force converge to respective steady-state values, typically within the sound crossing 
time over the distance equal to the orbital diameter or within about an orbital period when 
A4 is of order unity. The primary reason for this is of course because the perturbed density 
decreases quite rapidly with r and also because gravity is an inverse-square force. This is 
unlike the case of a straight-line trajectory where the drag increases secularly as ln[Vpt) for 
> 1 . The fast convergence of the drag force guarantees that one can use the steady-state 
values of X<^ and Tr for all practical purposes. 

Next, we check the dependence of the drag force on rinin(<^ Rp)- Figure 7 shows the 
results for Ai = 0.5, 2.0, and 4.0. The sizes of errorbars associated with finite grid resolution 
are smaller than the radius of a solid circle at each data point. First of all, the drag force, 
Tr, in the radial direction converges to a constant value as rmin decreases for both subsonic 
and supersonic cases. On the other hand, the drag force, X^, in the opposite direction 
of the orbital motion, varies as ln(l/rmin) with decreasing rmin for supersonic perturbers, 
while independent of rmin for subsonic perturbers. This dependence of on small rmin for 
circular-orbit perturbers is exactly the same as in the linear trajectory cases, which makes 
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sense because the curvature of a circular orbit is almost negligible in a tiny volume near the 
perturber. 

We plot in Figure 8 the steady-state drag forces for a circular-orbit perturber as functions 
of the Mach number. For all the points, Tmin = Rp/ 10 is taken and numerical convergence is 
checked. Filled circles give while open circles are for Xr/M.'^, which can be compared 

with Figure 3 of 099. For practical use, we fit the data using 

( 10 3-51A1-4.22^ < ^ -^ 

Xr^\ 0.5 In [9.33M2(_A^2 _ q 95^]^ f^j. 1.1 < < 4.4^ (13) 
i 0.3 for4.4<Al, 

and 

\ 0-7706 In (roooOTssAl) " 1-4703M, for M < 1.0, 

X^^l In [330( V^min)(M - 0.71)5-^2M-9-58] , for 1.0 < M < 4.4, (14) 
[ In [(i?p/r^in)/(0.11At + 1.65)] , for 4.4 < M, 

which are drawn as solid lines in Figure 8. The fits are within 4% of our semi-analytic 
results for Ai < 4.4 and within 16% for M. > 4.4. Both components of the drag force peak 
at ~ 1.2 — 1.4, analogous to the linear-trajectory cases, and X^p dominates over Xr for 
transonic perturbers. Although Xr is larger than X^ for Ai ^ 2.2, its effect on orbital decay 
of a perturber is insignificant (see §4). The local bumps in the drag forces at ~ 4.6 and 
7.8 are due to the self-overlapping of a wake tail explained in §3.1.1. 

Figure 8 also plots as dotted curves the results of 099 for drag force on hnear-trajectory 
perturbers. Despite the difference in the shape of the orbits, the agreement between X^ 
and Ostriker's formula is excellent for the subsonic case. This presumably reflects the fact 
that other than bending, the wake structure created by a subsonic circular-orbit perturber 
is not signiflcantly different from the linear-trajectory counterpart (see §3.1.2). Note also 
that Ostriker's formula for the supersonic case, for which Vpt = 2Rp is adopted to fit our 
numerical results, is overall in good agreement with X^ for a range of A^. It is remarkable 
that the subsonic and supersonic drag formulae (with Vpt chosen adequately in supersonic 
cases) obtained from perturbers moving straight still give a reasonably good estimate for the 
drag even on circular-orbit perturbers. 

4. SUMMARY AND DISCUSSION 

We have taken a semi-analytic approach to study the gravitational wake and the asso- 
ciated drag force on a perturber moving on a circular orbit in an infinite, uniform gaseous 
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medium. This work extends Ostriker (1999) who studied the cases with straight-hne trajec- 
tory perturbers. The circular orbit generally causes the wake to bend along the orbit and 
creates a trailing tail. For a subsonic perturber, the density wake has a weak tail and is dis- 
tributed quite smoothly (see Fig. 5). On the other hand, a supersonic perturber can catch up 
with its own wake, possibly multiple times depending on the Mach number, forming a very 
pronounced trailing tail across which density changes abruptly (see Fig. 3). Although the 
region influenced by the perturber keeps expanding with time into the surrounding medium, 
the drag force promptly saturates to a steady state value within less than the crossing time 
of sound waves over the distance equal to the orbital diameter. 

Because of asymmetry in the density wake, it is desirable to decompose the drag force 
into two components: and in the radial and azimuthal directions, respectively (eq. [11]). 
Figure 8 plots our main results for Xr and X^p as functions of the Mach number Al; equations 
(13) and (14) give the algebraic fits to the numerical results. The azimuthal drag force varying 
rather steeply with M. peaks at ~ 1.2 — 1.4, while the radial force is highly suppressed 
at A4 < 1 and exceeds the azimuthal drag aX }A ^ 2.2. It is remarkable that the drag on 
linear-trajectory perturbers given in 099 is almost identical to X^p for subsonic cases, and 
gives a good approximation for supersonic cases, too, provided Vpt = 2Rp. 

A striking feature in gravitational wakes formed by circular-orbit perturbers is a long 
tail in a trailing spiral shape. Such a tail structure is indeed commonly seen in recent 
hydrodynamic simulations for black hole mergers in a gaseous medium (e.g., Escala ct al. 
2004, 2005; Dotti ct al. 2006; sec also Sanchez- Salccdo & Brandenburg 2001). Albeit much 
weaker, it is also apparent in N-body simulations for satellite orbital decay in a coUisionless 
background (e.g., Weinberg 1989; Hernquist & Weinberg 1989). For supersonic perturbers, 
the tail is a curved Mach cone and bounded by the shock fronts that propagate radially 
outward. As explained in Appendices B and C, the outer edge of the tail is described by 
fit — ip = w+ + y+, where and are defined through equation (B2). By taking a time 
derivative of this equation for fixed (p, one can show that the propagation speed i?sh of the 
spiral tail in the radial direction is given by i?sh/cs = A1^(A^2^2_i)-i/2 for R = R/R^ > 1. 
Clearly, i?sh = for R/ Rp ^ 1, and R^h = CsA4/{A4^ — lY^"^ near R/ Rp = 1 as equation 
(CI) implies. This prediction is roughly consistent with Escala et al. (2004) who empirically 
found that the tail in a model with A4 = \/2 has an average propagation speed of Rsh ~ l.lc^. 

Many previous studies on dynamical friction adopted the drag formula based on per- 
turbers moving straight and estimated the Coulomb logarithm by taking rmax = -Rsys, where 
Rsys is the system size (see references in §1). Unlike a straight-line trajectory, a circular 
orbit naturally introduces a characteristic length scale, Rp. The results of our semi-analytic 
analyses suggest that, at least in a gaseous medium, the drag force obtained for linear- 
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trajectory perturbers can be a reasonable approximation to that for circular-orbit perturbers 
with the same Mach number if rmax = Vpt is taken equal to 2Rp. Since i?sys > 2i?p in 
most relevant situations, using rmax = -Rsys usually overestimates the drag force for ob- 
jects in near-circular motion. From hydrodynamic simulations for orbital decay of per- 
turbers in a stratified gaseous sphere, Sanchez-Salcedo & Brandenburg (2001) suggested 
ln(A) = ln(rmax/Tmin) —>■ 2.34 ln(0.79i?p/rmm) with an identification rmm = 2.25i?soft, where 
Rsoit is the softening radius of the point-mass potential (Sanchez-Salcedo & Brandenburg 
1999). This happens to be similar to our suggestion ln(A) — > ln(2i?p/rinin) for the parameter 
range Rp/rmin ~ 2 — 6 that they considered. 

Our suggestion for the Coulomb logarithm on near-circular orbit perturbers in a gaseous 
medium is also consistent with the common prescription for orbital decay of small satellites in 
coUisionless backgrounds (e.g., Tremaine 1976; Lin & Tremaine 1983; Hashimoto et al. 2003; 
Fujii et al. 2006). In particular, Hashimoto et al. (2003) found that Chandrasekhar's formula 
with ln(A) — > ln[it!p/(1.4i?soft)] gives excellent fits to the results of their N-body simulations. 
The discrepancy between our suggestion and their prescription may be attributable in part to 
the effects of coUisionless nature, nonuniform background density, and self-gravity of particles 
in their models. 

Density inhomogeneity in the background can also make significant changes to the clas- 
sical drag force. For coUisionless backgrounds, Just & Peiiarrubia (2005) found that nonuni- 
form density induces an additional drag force in the lateral direction of the pcrturbcr motion, 
taking up to 10% of the drag in the backward direction. They proposed the inverse of a 
local density gradient as an appropriate choice for rmax in the Coulomb logarithm (see also 
Spinnato et al. 2003), which was confirmed by Arena & Bertin (2007) who ran a number 
of numerical simulations for systems with a large central density concentration. Hydrody- 
namic models in Sanchez-Salcedo & Brandenburg (2001) studied the combined effects of 
nonuniform backgrounds and curvilinear orbits on dynamical friction in a gaseous medium, 
although it is challenging to isolate each effect separately. 

Finally, as a heuristic example, we consider the dynamical friction of a galaxy on a 
near-circular orbit subject to the drag force given by equations (13) and (14) due to an 
intracluster medium. The cluster is dominated by dark matter whose mass distribution 
follows the NFW profile with the characteristic mass Mo = 6.6 X 10 ^"^Mq and the scale 
radius Vg = 460 kpc (Navarro et al. 1996). The intracluster medium has a constant electron 
density rif. = 0.05 cm^^ and temperature 1 kcV; the corresponding adiabatic speed of sound 
is Cg = 500 km s~^. Initially, a galaxy with size rmm = 10 kpc and mass Mp = 5 x IO^^Mq 
including a dark halo (e.g., Zentner & Bullock 2003; Halkola et al. 2007) is orbiting at 
Ro — 100 kpc with an equihbrium velocity Vp — 10^ km s~^ about the cluster center. The 
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equations of motion in the orbital plane are 

R-R^^^ -{T/M,)In - (15a) 

R(p + 2Rip = -{T/Mp)!^, (15b) 

where $nfw is the NFW potential and dots represent derivatives with respect to time. 
Figure 9 plots the resulting temporal variations of the orbital radius of the galaxy. The solid 
line corresponds to the case with full 1^ and given by equations (13) and (14), while the 
dashed line is for a controlled case where is artificially taken to be zero. Except for slight 
mismatches in phase, both agree fairly well with each other, demonstrating that the radial 
drag force does not have a serious impact on the orbital decay. As equation (15) implies, 
it is the azimuthal drag that takes away most of the orbital angular momentum from 
the galaxy; the radial drag changes the orbital eccentricity more than the semi-major axis 
(Burns 1976; see also Just & Penarrubia 2005). 

Figure 9 also plots as dot-dashed hne the decay of the orbital radius blindly using 
Ostriker's formula with ln(T^i/rinin) = 4.6, corresponding to Vpt — 1 Mpc. While the galaxy 
motion remains supersonic {t < 1.5 Gyr), this value of the Coulomb logarithm overestimates 
the realistic drag force, on average, by a factor of 1.7 and thus brings the galaxy to R = O.li^o 
in 2 Gyrs, which is about 2.3 times faster than the estimate based on the realistic drag force. 
On the other hand, the result of using Vpt = 2R{t) in the formula of 099, shown as dotted 
line in Figure 9, is in excellent agreement with those using equations (13) and (14). This 
demonstrates again that Ostriker's formula with Vpt — 2R gives quite a good approximation 
to the drag force even in the case of a circular-orbit perturber. 

We are grateful to E. Ostriker for suggesting this project and making stimulating com- 
ments. We also acknowledge a thoughtful report from the referee, A. Escala, and helpful 
comments from J. Sanchez-Salcedo. This work was supported by Korea Science and Engi- 
neering Foundation (KOSEF) grant ROl-2004-000- 10490-0 at SNU. H. Kim has been partially 
supported by the BK21 project of the Korean Government. The numerical computations 
presented in this work were performed on the Linux cluster at KASl (Korea Astronomy and 
Space Science Institute) built with funding from KASI and ARCSEC. 



A. Solutions for V{x,t) Near the Perturber 

In this Appendix, we explore the approximate solutions to equation (8) near the per- 
turber (1^— 1|, \z\ ^ 1). For subsonic motion {M < 1), \wi\ ^ 1 and d 
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with 5 = R — 1. In this case, equation (9) has a single root, 

Inserting Wi into equation (8) and retaining the first-order terms in — 1|, and \z\, one 
obtains T>(x) = [s^ + (1 — A4'^){5'^ + at a steady state, whose functional form is 

identical to equation (9) of 099 for < 1. Obviously, the equi-density surfaces of this 
distribution are oblate spheroids with the short axis parallel to the direction of the motion 
of the perturber. 

For supersonic motion {Ai > 1), we seek for asymptotic solutions of V just behind the 
perturber along its orbit (i.e., < —s <^ 1, R = 1, z = 0). We further restrict the discussion 
to not-so-highly supersonic cases^, for which equation (9) has three roots: 

^ (A2a) 

/ A^^sinwo\~^ , 

Ws Wq — s I 1 , ( A2c) 

V Wo J 

as long as < — s <^ {M. — 1)^/^. Here, wq is a non-zero solution of 

M(2-2coswo)^/^ + Wo = (A3) 

for given M. I). Note that Wi and W2 are analogous to two solutions of 099 that contribute 
the wake in the rear Mach cone for a linear-trajectory perturber, while is a new solution 
introduced by a circular orbit. For M l"*"", wq —[2A{M — 1)]^^^- 

Substituting equations (A2) into equation (8) and taking the steady-state limit (7i — > 1), 
we obtain 

\s\ M^smwQ — Wo 

for —s<^ (Ai — 1)^/^. The first term in the right-hand side of equation (A4) again represents 
the perturbed density in the linear trajectory case. The second term, albeit small compared 
with the first term for small \s\, arises as the perturber on a circular orbit can overtake its 
wake. The effects of circular orbits are significant in regions away from the perturber, which 
is pursued numerically in §3. 



•^More precisely, 1 < M < Mi, where A^i is a critical Mach number discussed in Appendix B. 
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B. Number of Roots for Wi of Equation (9) 

To obtain the perturbed density V{x,t) at given time and location, it is crucial to find 
Wi that satisfies equation (9). Let us define f{w) = Aid{w; R, z) for a fixed set of Ai,R, z, 
and s. Then, finding solutions of equation (9) is equivalent to finding the intersections of 
two curves y — —w — s and y = f{w) on the w - y plane. Since \df /dw\ < M. for all R and 
it is trivial to show that there is only one solution for all space and time \l M. <1. This 
implies that any location inside the casual region for sound waves is influenced just once by 
a perturber moving at a subsonic speed. The same is true for a subsonic perturber moving 
on a straight-line trajectory studied by 099. 

A necessary condition for having more than one root is \df /dw\ > 1 for any w, which 
after some algebra results in 

M>1, R> M-\ and \z\ < M-^ {M^ - 1)'^^ {WR? - lf'\ (Bl) 

When the conditions (Bl) are fulflUed, there are two tangent points of the curve y — f{w) 
to lines with a slope of —1 in the range — tt < w < 0. Let {w+, y+) and {w-, y_) denote 
the tangent points, such that — tt < < w+ < and < < y_ < Al. The tangent 
condition df /dw\w± — —1 yields 

= -cos ^ J, (B2a) 

y± = \{M^R^-iy/^T{M^-lf^\, (B2b) 

at the z — plane. Note that the 27r periodicity of y = /('f^) ensures that {w± + 27m, y±) 
for arbitrary integer n are also tangent points to a line with slope of —1. 

For given M. > I, the number of solutions of equation (9) depends on the slopes of 
lines connecting the points {w^ — 2?™, y_) and {w+,y^). When the line passing through 
{w- — 27r,y_) and {w+,y+) has a slope larger than —1, which occurs for 1 < Al < A4i 
(see below for the definition of A4n with integer n > 1), the spatial regions that satisfy 
equations (Bl) as well as + y+ < —s < + have three roots for Wi, while the other 
regions have only a single root. When Aii < Ai < the slope of the line connecting 

(w_ —AT[,y_) and (w+,i/_|_) is greater than —1, so that the regions bounded by the condition 
Wj^ + |/+ < — s < W- + y- — 2tt possess five roots; the other regions have three roots. Note 
that the volume that does not satisfy conditions (Bl) still has a single root. As the Mach 
number increases further, new regions emerge to have a larger (odd) number of roots for 
Wi- Generalizing the discussion given above, the critical Mach number below which 
a steady-state wake possesses up to 2n -I- 1 roots can be determined by the requirement 



-17- 



(y+ — y-)/i'w+ — W- + 2?™) = —1. Using equations (B2), this is simplified to 

2{Ml - 1)^/2 - cos-i {2M-^ - 1) = 27in, (B3) 

where we take R — 1 without any loss of generality since the existence of multiple roots 
requires R > (eq. [Bl]). The first few solutions of equation (B3) are A4i — 4.6033, 

M2 = 7.7897, Ms = 10.9499, etc. 

The presence of multiple roots implies that the density wake in those regions are con- 
structed by a superposition of sonic signals that were emitted by the perturber at as many 
different locations (and corresponding retarded times) as the number of roots. The Ai- 
dependency of the number of roots for a circular-orbit supersonic perturber is in sharp 
contrast to the case of a straight-line trajectory where only one (within a sonic sphere) or 
two (within a rear Mach cone) roots are allowed regardless of the Mach number (099). When 
a perturber is moving at a supersonic speed along a circular orbit, it is able to catch up with 
its own wake (possibly multiple times), adding new sonic disturbances to the regions that 
were already disturbed by the perturber. Regions of multiple roots usually take a form of a 
long trailing spiral, as exemplified in Figures 2, 3, and 4. 



C. HALF-OPENING ANGLE 

Although a Mach cone is generally curved for a circular-orbit perturber, its half-opening 
angle in the vicinity of the perturber is the same as in the linear trajectory counterpart. As 
explained in §3.1.1, the interior of a curved Mach cone and an attached trailing tail is the 
region where equation (9) has multiple roots in a steady state. It is bounded by —s — w±+y± 
curves in the z — plane. More specifically, the —s — w+ + y+ curve describes the inner 
and outer edges of the Mach cone head near the perturber as well as the outer edge of the 
tail, represented by light solid curve and dashed curve at the boundaries of the shade region 
displayed in Figure 1. On the other hand, the —s = +y- curve corresponds to the inner 
edge of the tail (dot-dashed curve in Fig. 1) that meets the Mach cone head at R = Ai^^ . 
By Taylor expanding w+ and about i? = 1, s = in the z = plane and adding the 
resulting expressions together, one can show the head of the Mach cone is described by 

s ^ -{M' - iy/'\5\, (CI) 

up to the first order in \6\ = |i? — 1|. The negative sign in equation (CI) results since s is 
measured in the counterclockwise direction from the perturber (Fig. 1). Let 9 be the half- 
opening angle of the Mach cone near the perturber. Then, tan^ = = (Al^ — 1)"^''^, 
or ^ = sin~^(l/A^), the same result as in the linear trajectory case. 
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Fig. 1. — Schematic diagram illustrating the situation on the z — plane at time t. A 
perturber initially introduced at {Rp, 0, 0) moves along a circle with radius Rp at uniform 
angular speed in the counterclockwise direction, and is currently at the position Xp = 

(Rp, Qt, 0). At this time, an observer sitting at x = {R, ip, z) receives a sonic signal that was 
emitted by the perturber when it was at x' = {Rp, Qt', 0), where f = t — \x — x'\/cs is the 
retarded time. The angular variables arc w = ip' — Lp and s = Lp — Qt along the orbital plane. 
The shaded area represents a curved Mach cone and a wake tail formed by a supersonic 
perturber with Ai > 1. The inner edge (dashed curve) of the Mach cone meets the inner 
edge (dot-dashed curve) of the tail at a point on a circle with radius Rp/M.. The outer edge 
(hght sohd curve) of the Mach cone defines the outer edge of the tail. 
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Fig. 2. — Temporal evolution at the z = plane of (a) the dimensionless density wake 
D{x,t) in logarithmic color scale and (b) the number of roots for Wi in equation (9) for the 
case of A4 = 2.0. The time is in unit of Rp/cg. See text for details. 
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Fig. 3. — Density distributions of the steady-state wake for Ai = 2.0 on the x/Rp = 1 
{bottom right), y = {top left), and z = {bottom left) planes. The perturber is located 
at {x, y) = {Rp, 0), and the black circle in the bottom left frame denotes the orbit of the 
perturber. Color bar labels log P. 
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(a) M = 4,0 (b) M = 5.0 




Fig. 4. — Distributions of the perturbed density V in logarithmic color scale (top) and 
the corresponding number of roots for Wi in equation (9) (bottom) at the z = plane of 
the steady-state wake for (a) M. = 4.0 and (b) A4 = 5.0. The perturber is located at 
(x, y) = (Rp, 0). 
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Fig. 6. — Time evolution of the drag force for Ai = 0.5, 1.2, 2.0, 4.0, and 5.0. Solid curves 
draw the azimuthal drag X<^, while dashed curves are for the radial drag 1r. For all cases, 
'"min = Rp/ 10 is taken. Note that both and Ir converge typically within the sound 
crossing time across 2Rp, indicating that a steady state is reached quite rapidly. 
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Fig. 7. — Dependence on rmin of the steady-state drag force for M. = 0.5, 2.0, and 4.0. Left: 
For Rp/r^in > 10, the radial drag force 1r converges to a respective constant value marked 
by dotted line. Right: The azimuthal drag force varies as In (Rp/r^in) for small rmin when 
A4 > 1, while independent of rmin for A1 < 1. Dotted hnes indicate a slope of unity. 
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Fig. 8. — Gravitational drag force on a circular-orbit perturber in a gaseous medium as a 
function of the Mach number M.. The open and filled circles give the results of our semi- 
analytic calculation for the drag in the radial and azimuthal directions, respectively. For all 
the points, r^\^/Rp = 0.1 is taken. Solid lines draw the fits, equations (13) and (14), to the 
semi-analytic results. Dotted line corresponding to the force formula with Vpt = 2Rp in 099 
for the case of linear-trajectory perturbers is in quite a good agreement with the azimuthal 
drag for circular-orbit perturbers. 
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Fig. 9. — Orbital decay of a model galaxy caused by dynamical friction due to an intracluster 
gas. Solid line corresponds to the case when both Xr (eq. [13]) and (eq. [14]) are consid- 
ered, while dashed line shows the result with only 1^p (i.e., with Xr = 0). For comparison, 
the results of ©striker's formula with fixed Vpt = 1 Mpc and varying Vpt = 2R{t) are plotted 
as dot-dashed and dotted lines, respectively. 



